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I We consider the neutron matter quark matter phase transition along with possible 

existence of hybrid quark stars. The equation of state for neutron matter is obtained 

m 

using a nonperturbative method with pion dressing of the neutron matter and an 
analysis similar to that of symmetric nuclear matter. The quark matter sector is 
^ \ treated perturbatively in the small distance domain. For bag constant S 1 / 4 =148 

MeV, a first order phase transition is seen. In the context of neutron quark hybrid 
stars, Tolman-Oppenheimer-Volkoff equations are solved using the equations of state 
for quark matter and for neutron matter with a phase transition as noted earlier. 
Stable solutions for such stars are obtained with the Chandrasekhar limit as 1.58 
Mq and radius around 10 km. The bulk of the star is quark matter with a thin crust 
of neutron matter of less than a kilometer. 
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I. INTRODUCTION 



It is widely believed that neutron matter undergoes a phase transition to quark matter 
at high densities and/or high temperatures. The high temperature limit is expected to 
have interesting consequences in heavy ion collision and/or cosmology, whereas high baryon 
density behaviour is important for the study of neutron stars. 

The problem here for the treatment of nuclear matter or neutron matter equation of 
state is basically nonperturbative. A treatment of the same was developed by Walecka |IJ 
which consists of interactions of nucleons with a neutral scalar field o as well as a neutral 
vector meson to. These calculations however use meson fields as classical, and, use a a- 
field which is not observed. An alternative model for infinite nuclear matter consisting 
of interacting nucleons and pions was considered [@,|J recently where scalar isoscalar pion 
condensates simulate |4],|| the effects of a mesons with a phenomenological short distance 
repulsion arising from composite structure of nucleons or through vector meson exchanges. 
We extend this nonperturbative technique to neutron matter with the same parameters as 
for symmetric nuclear matter ||, and use this for hybrid stars. 

We first consider neutron matter equation of state at finite temperature through ther- 
mofield dynamics || as earlier || . A first order phase transition between neutron matter and 
quark matter seems to be indicated. As stated, neutron matter is treated in a nonpertur- 
bative manner. Quark matter is treated perturbatively for high densities or short distances 
EH- Solutions of the Tolman-Oppenheimer-Volkoff (TOV) equations yield, with the above 
equations of state, a hybrid star having a quark core and a crust of neutron matter, the 
two states of matter being in equilibrium at the interface ||. We obtain the Chandrasekhar 
limit for such hybrid stars as 1.58 M Q 

The outline of the paper is as follows. In Sec. II we shall consider neutron matter 
with pion condensates at zero temperature as well as finite temperature. In Sec. Ill we shall 
consider a standard form for quark matter with pressure given in terms of chemical potential 
and discuss about the phase transition from neutron matter to quark matter. Section IV 
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will consist of calculations of the properties of hybrid stars such as total mass, radius and 
moment of inertia. In section V we summarise the results. 



II. NEUTRON MATTER EQUATION OF STATE 



We shall consider here neutron matter equation of state at finite temperatures similar 
to the treatment of the same for symmetric nuclear matter ||. We start with the effective 
Hamiltonian for pion nucleon interactions H|J as 



H N (x) = Hk'(x) + H int (x) 



(1) 



where the free nucleon part Tt^ is given as 



(2) 



and the effective pion nucleon interaction part is given as 

iG _ _ G 2 



H int {x) = ip\(x) 



a. p<p + 



(3) 



In the above, e x = (M 2 — \J 2 X ) 1 ^ 2 where M denotes the nucleon mass. Furthermore the free 
meson part of the Hamiltonian is given as 

1 v, 



2 + (V0i) ■ (V&) + ™ 2 tf 



(4) 



Clearly in the above tpi refers to the large component of the nucleon spinor and <fi = Ti<pi. 

We expand the field operator <pi(z) in terms of the creation and annihilation operators 
of off-mass shell mesons satisfying equal time algebra as 

1 



(aj(F) 1 ' + a,i(z)) 



(5) 



and 



(6) 
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where we take with the perturbative basis cu z = (m 2 — S/l) 1 ^ 2 , with m denoting the mass of 
the meson. 

We shall now consider finite temperature neutron matter. For this purpose we shall use 
the methodology of thermofield dynamics |J. We shall have the temperature dependent 
background off-shell pion pair configuration given as 

\f,P>= Ui(0)\f >= Ui(P)U\vac > (7) 

where U is given as 

U = exp (- J f(k)<H{kyai(-k)Uk - h.c^j (8) 

and Ui(/3) describes the effect of temperature. The expression for this in terms of ordinary 
and thermal modes is given as |3|.[J 

U I (/3) = exp(B}(l3)-B J (P)) (9) 



with 



In the above, 



6i(fc)t = Uai{k)rf (11) 

so that bi(k)\f >= and 6s{k,l3) is a function to be determined later. The parallel unitary 
transformation as in Eq.(9) for the temperature dependance in neutron sector is given as |3j 

U II (P) = exp{B\ I {P)-B JI (p)) (12) 

with 

B n ((3y = ~J e F (k,P)MkyM-kydk (13) 
where 6p(k,/3) will be determined later. Thus we have the neutron matter density 



p{(3) = < vac\Uii(fiy^ a (xy^ a (x)Uii(l3)\vac > 

= 7(2tt)- 3 f dksin 2 6 F (14) 



where the degeneracy factor 7 = 2 for neutron matter ||. sin 2 0F is the distribution function 
for the fermions. Clearly, with sin 2 9p = Q(kf — k), the step function, Eq. (14) gives 
p = , ~fk'j/6TT 2 at zero temperature. For the interacting system sin 2 9p will be determined 
from the minimisation of the thermodynamic potential, as is done latter. 

We shall first calculate different contributions to the energy expectation values corre- 
sponding to the Hamiltonian as in Eq. (1) and (4). We thus have for the nucleon kinetic 
term 

h f {(3) = < vac\U n (pyM^^^M?)Uii((3)\vac > 

= 7 (2tt)- 3 J 'dk^szn 2 6 F (15) 

The kinetic energy due to the mesons is given by 

h k (J3) = <f,/3\H M (Z)\f,P> 

= 3(2vr)~ 3 y dkoo(k)[sinh 2 f(k)cosh2e B (k,(3) + stnh 2 6 B (k, (3)] (16) 



where ui{k) = y k 2 + m 2 . We next derive the interaction energy density from Eq. (3) as 

hint(P) = < f,mn((3yM^M^^ 2 ^)Uii((3)\f,/3> 

^2 



^p(P)<f,(3\: ( j> i (x)UZ)--\f,P> 
~2jT^ ( 17 ) 



where 



3 f dk sinh2f(k)cosh29B . ,97,^ , . ,9„ 1 

hu = 777^ / — -7 + sinh 2 f{k)cosh26 B + sinh 2 9 B . (18) 

{i^y J uj(k) y 2 1 

We shall now assume a phenomenological term corresponding to meson repulsion due to 
composite structure of mesons which is given as || 
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h*(P) = 3a(27r)- 3 / (sink 2 f{k)cosh29 B + sinh 2 9 B ) e^dk, (19) 

where a and are constants. Finally, the nucleon repulsion term is taken as || 

h R = \p 2 (P) (20) 

where p(/3) is as given in Eq. (14) and A is another constant. This repulsion term could 
arise from u;-exchange or otherwise || . 
Finally, the energy density is given as 

E((3) = (h f ((3) + h m (/3) + h R ((3))/p((3) (21) 

where 

h m W) = h k (f3) + h*(/3)+h int ((3). 
The thermodynamic potential density Q is given by [Hj 



n(P)=E((3)p-^-p BP (22) 

where the last term corresponds to nucleon number conservation with p B as the chemical 
potential. We may note that we shall be considering temperatures much below the nucleon 
mass so that in the expression for p(f3) we do not include antiparticle channel. The entropy 
density above is a = o B + o B with op being the entropy in fermion sector given as 

up = — ^ / dk sin 2 9p(k, j3)ln{sin 2 0p(k ) (3)) 



(2vr) ; 



+ cos 2 6 F (k, (3)ln{cos 2 6 F {k, (5)) 



3 

o B = — - — -r / dk sinh 2 8 B (k, f3)ln(sinh 2 6 B (k, (3)) 



and similarly the meson sector contribution o B is given as 

3 

W) 

- cosh 2 6 B (k } (3)ln(cosh 2 6 B (k, (3)) 

Thus the thermodynamic potential density Q is now a functional of 9p(k,l3), 8 B (k,(3) as 
well as the pion dressing function f(k) which will of course depend upon temperature. 
Extremisation of Eq (21) with respect to f(k) yields 
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tanh2f(k) = % =, - 2 =, . (23) 

Similarly minimising the thermodynamic potential with respect to 0B(k,f3) we get 

sinhH B = (24) 

where 

J = {U + 2MJ + aeRlk2 ) cosh2 hk) + ^£sinh2f(k) (25) 
Once we substitute the optimised dressing as in Eq. (23), the above simplifies to 

J = r u + 9lP + ae R * k y/\u + ae^) 1 ' 2 (26) 
Moj 

which is different from uo due to interactions. Further, minimising the thermodynamic 
potential with respect to 9p{k,P) we have the solution 

sin 2 9 F = —, (27) 

with 

tF = Wf hM + 2 * p + &f (28) 

where 7 2 m is given in equation (18). We may note that the change in ep above from k 2 /2M 
is also due to interaction. 

To calculate different thermodynamic quantities as functions of baryon number density 
we first use Eq. (14) to calculate the chemical potential in a self consistent manner with 
p and pb occuring also inside the integrals through sin 2 9p as in Eq. (27). Thus for each 
p, we determine ps so that Eq. (14) is satisfied. The ansatz functions 9b, 9f and f(k) get 
determined analytically through the extremisation of the thermodynamic potential. The 
parameters a, A and R n are fixed so as to reproduce the ground state properties of zero 
temperature nuclear matter as is done earlier ||. The parameter values A = 0.54 fm 2 , 
R n = 1.18 fm and a = 0.12 GeV give the nuclear matter single particle energy as —15.03 



MeV at the saturation density of 0.19 fm~ 3 . We shall also use here the same values for a, 
and A. 

With the thermodynamic potential determined as above, we calculate different thermo- 
dynamic quantities. Pressure is calculated using the thermodynamics relation |T(J 

P((3) = -n(0) (29) 

Pressure as a function of neutron matter density is plotted in Fig. 1 at temperatures of 0, 
10, 15 and 20 MeV. We note that as earlier [|2]|| the equation of state is quite soft. In Fig 2 
we plot pressure versus baryon chemical potential at temperatures of 0, 5 and 10 MeV. We 
note from Fig. 1 that at zero temperature pressure for neutron matter vanishes at number 
densitiy of about 0.1/fm 3 so that the pressure at the surface for a neutron star can be zero 
with a finite number density. This unusual feature is due to the presence of nonperturbative 
contributions in equation (28). 

In the next section we shall consider quark matter as well as the transition from neutron 
phase to quark phase. The zero temperature version of equation (29) shall be used for phase 
transition studies. 



III. QUARK MATTER EQUATION OF STATE AND PHASE TRANSITION 

Existence of quark matter in the core of neutron stars/pulsars is an exciting possibility 



TTHl3|| . Densities of these stars are expected to be high enough to force the hadron con- 
stituents or nucleons to overlap thereby yielding quark matter. Since the distance involved 
is small, perturbative QCD is used to derive quark matter equation of state. We take the 
quark matter equation of state as in Ref JFj in which u,d and s quark degrees of freedom 
are included in addition to electrons. As here [0] we set the electron, up and down quark 
masses to zero and the strange quark mass is taken to be 180 MeV. In chemical equilibrium 
fJ>d — fJ's = A*« + A*e- In terms of baryon and electric charge chemical potentials fis and [Ae, 
one has 
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1 



1 1 



1 



1 



(30) 



The pressure contributed by the quarks is computed to order a s = g 2 /^ where g is 
the QCD coupling constant. Confinement is simulated by a bag constant B. The electron 
pressure is 

,4 



e 12vr 2 ' 

The pressure for quark flavor f, with f=u,d or s is 

1 



(31) 
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4tt 2 

7T 3 



- 2.5m}) + 1.5mJZn ' M/ + kf 



fjLfkf — mj/ri 



'n f + kj 



nit 



k 4 



(32) 



The Fermi momentum is kf = (fjjj — m 2 ) 1 ^ 2 . The total pressure, including the bag constant 
B, is 



P = P e + Y.Pf-B. 

f 



(33) 



There are only two independent chemical potentials \ib and \ie- He is adjusted so that the 
matter is electrically neutral, i.e. dP/dfiE — 0. The baryon number density is given by 
p = dP/d/j, B . 

We now consider the scenario of phase transition from cold neutron matter to quark 
matter. As usual, the phase boundary of the coexistence region between the neutron and 
quark phase is determined by the Gibbs criteria. The critical pressure and critical chemical 
potential are determined by the condition 



Pnmil^B) — Pqm\HB)- 



(34) 



The r.h.s of equation (34) is the same as P in equation (33), where as the l.h.s is the zero 
temperature limit of equation (29). We take a s = 0.4 and the bag constant B = (148 MeV) 4 , 
which is a reasonable value to calculate pressure in the quark sector. In Fig 3 we plot pressure 
versus chemical potential hb for cold neutron matter and quark matter as considered here. 
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The (P, fi) curves for quark matter and neutron matter yield the critical parameters as 
P cr = 4 MeV/fm 3 and (hb)ct — 960 MeV. The corresponding critical energy densities 
for neutron matter and quark matter respectively are given as e"™ = 226 MeV/fm 3 and 
ej?™ = 284 MeV/fm 3 . The corresponding baryon number densities are = 0.24/fm 3 
and pg 1 = 0.33/fm 3 . Thus there is a first order phase transition. We also note that the 
phase transition seems to occur around number density of only about one and half times the 
nuclear matter density. There is a second crossing point in (P, /i) curves not seen in Fig.3. 
However, it occurs at higher densities where we do not expect neutrons to exist because of 
overlap of quark wave functions so that neutron matter does not exist. 

The phase transition as above was possible due to the first two terms in equation (28). 
These terms arise in a nontrivial manner from interactions, with I2M given in equation in 
(18). We note that here we are using the same parameters for neutron matter as earlier for 
symmetric nuclear matter, giving rise to a soft equation of state ||. 

The early phase transition from neutron matter to quark matter obviously implies that 
the interior of "neutron star" will usually consists of quark matter. We investigate this 
possibility in the next section. 



IV. HYBRID STARS 

For the description of neutron star, which is highly concentrated matter so that the metric 
of space-time geometry is curved, one has to apply Einstein's general theory of relativity. 
The space-time geometry of a spherical neutron star described by a metric in Schwarzschild 
coordinates has the form fl4] , |15l 

ds 2 = -e u(r) dt 2 + [1 - 2M{r)/r]- 1 dr 2 + r 2 [d0 2 + sin 2 Qd<j) 2 } (35) 

The equations which determine the star structure and the geometry are, in dimensionless 
forms [ [IHI 



dP(fr ) _ ^ [e(fr ) + P(fr )] [M(f r ) + 47raf 3 P(fr )] 
df ~ ' f 2 [l-2GM{fr )/f] 
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rro 

M(fr ) = 4:iia I df'r' e(fV ), (36b) 



o 



and the metric function, z/(r) is given by 

dv(rr ) _ g [M(f r ) + 47raf 3 P(fr )] 
~~ r 2 [l -2GM(rr )/f] ' 

In equations (36) the following substitutions have been made. 

e = e/e c , P = P/e c , r = r/r , M = M/M Q , (37a) 
where, with /i = 197.329 MeV fm and r = 3 x 10 19 fm, we have 

a = e c r o 3 /M , G = (37 b) 



In the above, quantities with hats are dimentionless. G in equation (|37b|) denotes the 
gravitational constant with G = 6.707934 x 10" 45 MeV -2 . 



In order to construct a stellar model, one has to integrate equations ( |36a| ) to (|36c| ) from 



the star's center at r=0 with a given central energy density e c as input until the pressure 
P(r) at the surface vanishes. As stated in the last section, with any reasonable central 
density, we expect that at the center we shall have quark matter, and not neutron matter. 
Hence we shall be using here the equation of state for quark matter through equations (32) 
and (33) with P(0) = P(e c ). We then integrate the TOV equations until the pressure and 
density decrease to their critical values, so that there is a first order phase transition from 
quark matter to neutron matter at radius r = r c . For r > r c , we shall have equation of state 
for neutron matter where pressure will change continuously but the energy density will have 
a discontinuity at r = r c . TOV equations with equation of state for neutron matter shall be 
continued until we reach the density around 0.1 particles/fm 3 when the pressure vanishes 
(see Fig 1). This will complete the calculations for stellar model for hybrid "neutron" star, 
whose mass and radius can be calculated for different central densities. 

We now plot in Fig 4 the mass of a star as a function of central energy density to examine 
the stability of such stars. As may be seen from the figure dM/de c becomes negative around 

11 



1540 MeV/fm 3 after which they become unstable and may collapse into black holes with 



the Chandrasekhar limit as 1.58M JT^JT^]. The dashed part of the curve indicates the 
instability region. Fig 5 shows the mass as a function of radius for such stars obtained for 
different central densities varying in the range 900 to 1500 MeV/fm 3 . This yields stable 
hybrid stars of masses M ~ 1.49 to 1.58 M Q with raddi R ~ 9.85 to 9.4 km respectively, 
similar to the results of Ref. ||13|| . 



The energy density profile obtained from (|36a|) to (|36c| ) are plotted in Fig 6 for central 
densities e c = 900 MeV/fm 3 and e c = 1500 MeV/fm 3 . As we go away from the core towards 
the surface through TOV equations, when the critical pressure is reached, the density drops 
discontinuously indicating the first order phase transition. Thus e.g. for central density of 
900 MeV/fm 3 such a star has a quark matter core of radius 9 km with a crust of neutron 
matter of about 0.85 km, whereas for e c =1500 MeV/fm 3 , the quark matter core radius is 
8.9 km with neutron matter crust of 0.5 km. 

It may be interesting to consider smaller central densities, when neutron matter is ex- 
pected to be more abundant. For example, we have noted that at central densities of e c =700 
Mev/fm 3 to 350 Mev/fm 3 , R changes from 9.87 km to 7.3 km, whereas the quark core 
changes from 9 km to 5.1 km. The neutron crust here increases from 0.87 km to 2.2 km, 
and the mass of the star decreases to as low as 0.32M Q . It is not clear to us whether, such 
small hybrid stars or for that matter neutron stars can get formed. 

We may also calculate the surface gravitational red shift Z s of photons to see the possi- 
bility of distinguishing such a hybrid star from a neutron star or a quark star. It is given by 



1. (38) 



[1 - 2GM/R) 

In Fig 7 we plot Z s as a function of M/M & . As may be seen, Z s as related to mass and 
radii of hybrid stars is similar to that of pure neutron stars and hence a measurement of the 
same will not be an evidence for the existance of such a star. It is however possible that the 
discrete slowing down of pulsars due to the presence of two states of matter with various 
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radii shall throw some light on the above structure. For a quantitative estimate of the same, 
however, we shall need interactions at the interface, which at present looks impossible to 
tackle. This will need for example the expression for the moment of inertia of pulsars as 



well as the equilibrium for the rotating hybrid stars given as [18 



3 Jo y/[l-2GM(rr )/r] 
We may further take the relativistic Keplerian angular velocity Qk to be given as [T9j 



tt k M/M & , . 

0.72. 1 ° (40) 



lO^- 1 ' \ (R/lOkmf 
as for neutron stars and estimate the same. Fig 8 shows the variation of Qk as a function of 
M/Mq for such a star. 

In Fig 9 we plot the fraction of the quark matter in a star as a function of the central 
density. As can be seen from the figure quark matter constitutes the dominant contribution 
to the mass of such stars. The internal quark matter structure of such hybrid stars could 
influence neutrino emissivity due to larger number of /3-decay channels in quark matter 
phase which may be relevant for the dynamics of supernova explosions giving rise to pulsars, 
now expected to be really hybrid stars. 



V. CONCLUSIONS 

Let us summarise the findings of the present paper. 

Using a nonperturbative equation of state for neutron matter and a usual perturbative 
equation of state for quark matter we saw that a first order phase transition exists between 
the neutron phase and quark phase. We note that neutron matter equation of state practi- 
cally has no free parameters since the parameters are fixed || to yield the standard nuclear 
matter properties. Quark matter equation of state has the parameters : m s , a s and B. The 
bag constant B, turns out to be a very sensitive parameter. We found, for example, for 
a s = 0.4, neutron matter quark matter transition occurs for B ~ (148MeV) A . As can be 
seen in Fig. 3 there is not much freedom for B. 
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The phase transition from neutron matter to quark matter seems to indicate that the 
core of a "neutron star" shall consist of quark matter. To study the hybrid star, namely 
a star consisting of both quark matter and neutron matter, we applied the TOV equation 
to the appropriate equations of state with a given central energy density e c . It turns out 
that stable hybrid stars with a core of quark matter and a crust of neutron matter can exist 
upto e c ~ 1500 MeV/fm 3 beyond which instability is indicated. The bulk of the hybrid 
star is provided by the quark matter, the neutron matter providing only a thin crust. It is 
thus the quark matter equation of state and the corresponding parameters which shall play 
a dominant role in the formation of hybrid stars. Pulsars are expected to be stars of this 
type, but the gross properties appear to be similar to what we believe regarding neutron 
stars. We may recall that the nuclear matter equation of state used here is rather soft with 
compressibility of about 150 MeV. 

We have also calculated other gross properties of such hybrid stars like moment of iner- 
tia, surface gravitational red shift and Keplerian rotation period. The results are similar to 
those obtained by others [13]. 
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FIGURES 

FIG. 1. Pressure P as a function of neutron matter density p is plotted for temperatures T=0 
MeV, 10 MeV, 15 MeV and 20 MeV. 

FIG. 2. Pressure P of neutron matter as a function of baryon chemical potential is plotted for 
temperatures T=0 MeV, 5 MeV and 10 MeV. 

FIG. 3. We plot here pressure versus baryon chemical potential for neutron matter and for 
quark matter with a s =0.4 and B 1 ^ =148 MeV at zero temperature. The critical parameters may 
be seen as P cr =4 MeV/fm 3 and (fi B )cr=9Q0 MeV. 

FIG. 4. We plot mass of the hybrid star as a function of central density. The instability region 
is shown as dashed curve. 

FIG. 5. We plot mass as a function of radius for hybrid stars here. The dashed line here 
indicates the instability of the solutions. 

FIG. 6. We plot here energy density profile for central densities e c =900 MeV/fm 3 and e c =1500 
MeV/fm 3 . The discontinuities at the ends indicate a first order phase transition at the interface of 
quark matter and neutron matter. 

FIG. 7. We plot here surface gravitational red shift Z s versus the star mass. 

FIG. 8. We plot here Keplerian angular velocity Ofc versus star mass. 

FIG. 9. We plot here the ratio between the mass in the quark matter phase and the total mass 
of the star as a function of its central density. The dashed curve indicates the instability region. 
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